arXiv:1508.04767vl [math.NA] 19 Aug 2015 


On a family of Weierstrass-type root-finding methods with 

accelerated convergence 

Petko D. Proinov, Maria T. Vasileva 

Faculty of Mathematics and Informatics, University of Plovdiv, Plovdiv 4000, Bulgaria 


Abstract 

Kyurkchiev and Andreev (1985) constructed an infinite sequence of Weierstrass-type iterative 
methods for approximating all zeros of a polynomial simultaneously. The first member of this 
sequence of iterative methods is the famous method of Weierstrass (1891) and the second one is 
the method of Nourein (1977). For a given integer N > 1, the Mh method of this family has the 
order of convergence N + 1. Currently in the literature, there are only local convergence results 
for these methods. The main purpose of this paper is to present semilocal convergence results for 
the Weierstrass-type methods under computationally verifiable initial conditions and with compu¬ 
tationally verifiable a posteriori error estimates. 
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1. Introduction and preliminaries 

Throughout this paper (K, | • |) denotes an algebraically closed normed field and IK[z] denotes 
the ring of polynomials (in one variable) over IK. We endow the vector space IK" with the p-norm 
IWIp = (Z"=i \xi\ p ) 1/p for some 1 < p < oo, and we equip (M", || • || p ) with coordinate-wise ordering 
< defined by 

x <y if and only if x t < y t for each i - 1,..., n. (1.1) 

Then (M", || • || p ) is a solid vector space. Also we define a cone norm || • || in IK" with values in M" 
by 

Ifill = (fill, • • • , finl)- (1-2) 

Then (IK", H • ||) is a cone normed space over R" (see, e.g., Proinov [|9l). 

Let / 6 K[z] be a polynomial of degree n > 2. A vector if e IK" is called a root-vector of / if 
f(z) = a {) UUz - ifi) for all z e IK, where a 0 6 IK. 
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In 1891, Weierstrass l!T8l published his famous iterative method for simultaneous computation 
of all zeros of /. The Weierstrass method is defined by the following iteration 


Sk+l) 


x® - Wf{ x w ) 


ikh 


k- 0,1,2,.... 


(1.3) 


where the operator Wf : T> c IK" —» IK" is defined by 
W f (x) = (W l (x),...,W 1 (x)) with W t (x) = 


/(*«) 


ao Yljti (*« - *;) 


(f = 1, • • •, n), 


(1.4) 


where a 0 6 IK is the leading coefficient of / and the domain !) of Vfy is the set of all vectors in IK" 


with distinct components. The Weierstrass method ( |1.3[ ) has second-order of convergence provided 
that all zeros of / are simple. Other iterative methods for simultaneous finding polynomial zeros 
can be found in the books BHUHE, 15J and the references therein. 

In 1985, Kyurkchiev and Andreev f3]| introduced a sequence of iterative methods for approx¬ 
imating all zeros of a polynomial simultaneously. The first member of their family of iterative 
methods is the Weierstrass method ( |1.3[ ) and the second one is the method of Nourein [IT). 

Before we present Kyurkchiev and Andreev’s family of iterative methods, we give some nota¬ 
tions which will be used throughout the paper. We define the binary relation # on IK" by 


x#y <=> x, ±yj for all i,j e l n with i ± j. (1.5) 

Here and throughout this paper, we denote by I n the set of indices 1,...,/?. For two vectors x e IK" 
and y e R" we define in R" the vector 


X _ /jXlj _ _ _ \Xn\\ 

y \ yi ’ ’ y n ) 

provided that y has only nonzero components. We define the function d: K' ! R" by 

d(x) = (d i(x),..., d n (x)) with dj(x) = min \x t - x\ (i = 1,..., n). 

j*i 

In the sequel, for a given vector x in IK", x, always denotes the /th component of x. In particular, if 
F is a map with values in IK", then F,(x) denotes the /th component of the vector Fix). 

Definition 1.1. Suppose / 6 K[z] is a polynomial of degree n > 2. We define the sequence 
( T tN> )“ =( , of functions T <N} : D N c IK" —> IK" recursively by setting T (0> (x) = x and 

T\ n+ '\x) = x ( - f{Xi \ - (/ = 1, • ■ • ,n), (1.6) 

Y\(x, - T™(x)) 

j*i 

where the sequence of the domains D N is also defined recursively by setting Z) 0 = IK" and 

D N+l = {x e D n : x#T {N) (x)}. (1.7) 
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Let N 6 N be fixed. Then the Nth method of Kyurkchiev-Andreev’s family can be defined by 
the following fixed point iteration 


r (k+ 1 ) 


T (N \x {k} ), k = 0,1,2, 


( 1 . 8 ) 


Currently in the literature, there are only local convergence results for the Weierstrass-type 
methods Q (see (3l m. In this paper, we present semilocal convergence results for the 
Weierstrass-type methods under computationally verifiable initial conditions and with computa¬ 
tionally verifiable a posteriori error estimates. These results are obtained by using some results of 
mon and AD. 

The paper is structured as follows: In Section[2j we obtain new local convergence results (The¬ 
orem 2.11, Corollary 2.12 and Corollary |2.13[ ) for the Weierstrass-type methods ( |1.8[ ). In the case 
N = 1 (Weierstrass method) and p = oo the main result of this section reduces to a result of Proinov 
IfTOl Theorem 7.3]. In Section[3j we present our semilocal convergence results (Theorem |3.2[ The- 
orem |3.5[ Corollary |3.3| and Corollary |3.6[ ) for the Weierstrass-type methods ( |1.8[ ). Note that these 
results are based on the local convergence results obtained in the previous section. In Section [4] we 
provide three numerical examples to show the applicability of our semilocal convergence results. 

Throughout this paper, we follow the terminology of IfiOll . In particular we refer to this paper 
for the definition of the following notions: quasi-homogeneous function of degree r > 0; gauge 
function of order r > 1; function of initial conditions of a map; initial point of a map; iterated 
contraction at a point. 

2. Local convergence analysis of the Weierstrass-type methods 

Let / 6 K[z] be a polynomial of degree n > 2. In this section, we study the convergence of 
the Weierstrass-type methods ( |1.8[ ) with respect to the function of initial conditions E : T> —> M+ 
defined by 


E(x) 


x-f; 


d(x) 


for some 1 < p < oo. We define the function T: 

T(t) = (1 +2r)(l + 


by 


(n - 1 )p 


n— 1 


( 2 . 1 ) 


( 2 . 2 ) 


Throughout this section we denote by R the unique positive solution of the equation 'F(t) = 2. It 
can be proved that 

n(2 l/n - 1) 1 

(23) 

where q is the conjugate exponent of p, i.e. q is defined by means of 1 < q < oo and l/p + l/q = \. 
The lower estimate in ( |2.3[ ) can be proved as Lemma 7.4 of IflOl and the upper estimate is trivial. 
Also, we define the functions to: M + —» M + and cf>: [0,/?] —» [0,1] by 

to(t) - 1 


co(t) = 1 + 


(n - \)p 


2-1 


and (p{t) 


1 - 2 tco(t) 


(2.4) 




























It follows from the definition of 7? that 


m(7?) = 2/(1 + 2 R), oj(R) < 2, tu(7?) < 1/(2 R) and 0(7?) = 1. (2.5) 

Definition 2.1. We define the sequence (0,v)/:’ =o of nondecreasing functions <p N \ [0,7?] —» [0,1] 
recursively by setting 0 o (O = 1 and 


0V+1(0 


a>v(0 - 1 
1 — 2t(o N (t) 


where «»(<) = (l + 


( 2 . 6 ) 


Proof of the correctness of Definition 2.1. We prove the correctness of the definition by induc¬ 
tion. For AT = 0 it is obvious. Assume that for some N > 0 the function </> v is well-defined and 
nondecreasing on [0,7?] and cp N (R) = 1. We shall prove the same for (f> N+ 1 . From the induction 
hypothesis, we get that co N is nondecreasing on [0,7?] and a> N (R) = oMR). Then it follows from 
( [23) that 

1 - 2 toj N (t) > 1 - 2Rto N (R) = 1 - 2Rcjl>{R) > 0 (2.7) 


which guarantees that the function f N+1 is well-defined on [0,7?]. Obviously, <p N+x is nondecreasing 
on [0,7?]. From the definition of <p N + i, co N (R) = a>(R) and ([23), we get 

to N (R) — 1 ru(7?) — 1 


<pN+ l(R) 


m = i. 


1 - 2Rco n (R) 1 - 27?m(7?) 

This completes the induction and the proof of the correctness of Definition al [ 
Definition 2.2. Given N > 0, we define the function ip N : [0,7?] —> [0,7?] by 


□ 


<Av(0 = t(p N (t). 


( 2 . 8 ) 


Lemma 2.3. Let N > 0. Then: 

(i) (pN is a quasi-homogeneous of degree N on [0,7?]; 

(ii) (pN+ft) < (p{t)(p N {t)for every t e [0,7?]; 

(iii) 0v(O < <pif) N for every t 6 [0,7?]; 

(iv) ip N is a gauge function of order N + 1 on [0,7?]. 


Proof We prove Claim (i) by induction. The case TV = 0 is obvious. From the induction hypoth¬ 
esis and Example 2.2 of IfTUII , we conclude that the function oj N (t) - 1 is a quasi-homogeneous of 
degree TV + 1 on [0,7?]. Hence, f N+ \ is quasi-homogeneous of degree TV + 1 on [0,7?] as a product 
of a quasi-homogeneous of degree AT + 1 and a nondecreasing function. This ends the proof of (i). 
From the fact that the function oj^(t) - 1 is a quasi-homogeneous of degree TV + 1 on [0,7?], we 
obtain 

(o N (t) - 1 < f N (t)(<jj(t) - 1 ). 


From (|2.6[), the last inequality and co N (t) < of t), we get 


0V+1(0 


M - v(,) ~ 1 < M» " < ' ) " 1 

1 - 2 tojff(t) 1 - 2 tco{t) 


= (p N (t)(p{t) 


which proves (ii). Claim (iii) follows from (ii) by induction. Claim (iv) follows from (i) and 
definition of tp N . □ 
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Definition 2.4. For a given integer TV > 1, we define the increasing function : [0,7?] —> [0,1) 
by 

flh'(t) = <^v-i(0 - 1 (2.9) 

and we define the decreasing function if / N : [0,7?] —» (0,1] by 

<Av(0 = 1 - (2.10) 


where the function oj,y is defined in (2.6). 


Proof of the correctness of Definition 2.4. The functions fi N and if/ N are well defined on [0,7?] 
since ojn-i is well defined on this interval. The monotonicity of these functions is obvious. It 
remains to prove that fi N (R) < 1 and f N (R) > 0. It follows from a» w _i(7?) = co(R) and (|2.5[) that 


/3 n {R) = <u(7?) - 1 < 1 and fi N (R) = 1 - 2Rco(R) > 0 


which completes the proof of the correctness of Definition |2.1| fl 

Lemma 2.5. Let N > 1. Then 


(i) /3n is a quasi-homogeneous of degree N on [0,7?]; 

(ii) [3 N (t) = f N (t)f N (t) for every t 6 [0,7?]. 

Proof Claim (i) follows from Lemma |2.3[ i) and Example 2.2 of [ITOl . Claim (ii) follows from 
Definition 12.11 and Definition 12.41 □ 


Lemma 2.6 (lfT3l). Let f G K[z] be a polynomial of degree n > 2. Assume that ^ 6 IK" is a root- 
vector of f and N > 1. If x 6 D N+h then for every i e I„ 


T\ N+x \x)-fi 


i - y [(i+ u .j) 


j*i 


(Xi ~ fi). 


( 2 . 11 ) 


where uj e IK is defined by 


Uj = 


T ( f\x)-£j 
Xi - T { p(x) 


( 2 . 12 ) 


Lemma 2.7 (031). Let u,v,fe IK", a > 0 and 1 < p < oo. If y is a vector with distinct components 
such that 


\\u - f\\ < a||v-^l|, 


(2.13) 


then for all i, j e /„, 


Vj\ > 1 - (1 + or) 




d(v) 


(2.14) 


Lemma 2.8 (HU). Let u,v,f £ IK", a > 0 and 1 < p < oo. If y is a vector with distinct components 
such that ( |2.13[ ) holds, then for all i, j £ 


Uj - Uj\ > 


|l -2 1/? (1 +a) 



I v i ~ v j\ 


(2.15) 
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Lemma 2.9. Let f £ K[zJ be a polynomial of degree n > 2, f £ K" be a root-vector of f, N > 1 
and 1 < p < oo. Suppose x £ IK" is a vector with distinct components such that 


E{x) < R, (2.16) 

where the function E is defined by Then f has only simple zeros in IK, x 6 D N and 

||r w (x) - £|| </3 N (E(x))\\x-a (2.17) 

Proof It follows from Proposition 5.3 of [fTUl that the vector f has distinct components, which 

means that / has only simple zeros in IK. Further, we proceed by induction. If N = 1, then the 

proof can be found in IflOl . Assume that both a; 6 D N and ( ]2. 1 7[ ) hold for some N > 1. 


taking into account (2.16) and the fact 


First, we prove that a: 6 D N+l . It follows from ( ]2.17| ) that condition ( ]2. 1 3[ ) is satisfied with 
u = 7’ ( V, (a'), v = x and a = 1. Therefore, by Lemma 2.7 
that a; is a vector with distinct components, we obtain 


T[ n \x) | >11-2 


d{x) 


Xj\ > (1 - 2 E{x))dfx) > 0 


(2.18) 


for every j 4- i. Consequently, x#T iN \x) which proves that a: 6 D N+l . 

Second, we shall prove that ( |2.17[ ) is true for N + 1. Obviously, the last statement is equivalent 
to 

|T ; (A,+1, (a;) - fi\ < fi N+] (E(x))\xi - fi\ for every i e /„. (2.19) 

Let i £ I„ be fixed. We consider the vector u = (Uj) j± l £ IK" -1 , where Uj is defined by ( |2. 1 2[ ). It 


follows from (2.17) and (2.18) that 


\ T f\x)-U < fi N (E(x)) \xj-tj\ 


Tf\x)\ 


1 - 2 E{x) dj(x) 

Taking the p-norm and using Lemma [13} ii) and ffit) < 1, we get 

E(x)fi N (E(x)) E(x)f N (E(x))fi N (E(x)) 


< 


1 - 2 E{x) 1 - 2 E{x) 

Then, from Lemma [276} Lemma 3.2 of [IT2l . ( |2.20 [ ) and ( |2.9| ), we obtain 


< E{x)4> n {E{x)). 


( 2 . 20 ) 




]Td +«;)-! 


< 


1 + 


\Xi - fi\ 


n -1 


l 


I Xi - fi\ 


(n - \ f lp t 

(m N (E(x)) - 1)|a:/ - = p N+1 (E(x))\Xi - U 


which proves that (|2.19[) is true for N + 1. This completes the induction and the proof. 


□ 
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Lemma 2.10. Let f G K[z] be a polynomial of degree n > 2, £ e IK" be a root-vector of f, N > 1 
and 1 < p < oo. Let T (N) : D N c IK" —» IK" and E : D N — * M+ be defined by Definition 1.1 and 
( |2.1[ ), respectively. Then: 

(i) E is a function of initial conditions of T iN) with a gauge function <p,y of order N + 1 on 
J = [0, R\; 

(ii) T (N} is an iterated contraction at £ with respect to E with control function [i N ; 

(iii) Every point x ({)) 6 IK" such that E(x a>> ) G J is an initial point ofT iN) . 


Proof, (i) First, we prove that 

E{T (N \x)) < cp N (E(x)) for all a: 6 ID such that E(x) e J. 


( 2 . 21 ) 


Condition ( |2.17[ ) allow us to apply Lemma 28 with u = T (N \x), v = x and a = fi N (E(x)). There¬ 
fore, we get 

\T ( : N) (x) - T ( . n \x) 1 > (1 - 2 l/q E(x)(l +ME(xmxi - x,\ > f N (E(x))\ Xi - xA. 


Taking minimum over j t i on both sides of this inequality, we obtain 

dfT {N \x)) > f N (E(x))dfx) > 0. 

It follows from ( |2.19| ), ( |2.22[ ) and Lcmma [23| ii) that 

\*i ~ U 


( 2 . 22 ) 


\T\ N \x)~U ^(3 n (E(x)) \ Xi -&\ 

< -= <Pn(E(x))- 


dfTWfx)) ~ f N (E(x)) ddx) 

Taking the p-norm on both sides of this inequality, we get 

E(J iN \x)) < f N (E(x))E(x) = cp N (E(x)) 


d,(x) 


(2.23) 


which proves ( |2.21[ ). Now Claim (i) follows from ( |2.2 1 [ ) and Lemma [23| iv). 

(ii) It follows from Lemma [2791 

(iii) It follows from Lemma [29] that x i(y> e D N . According to Proposition 2.7 of [ITOil to prove 
that x i0> is an initial point of T (N) it is sufficient to prove that 


x G D/y and E(^x) G J => T^fx) G Dp . 


(2.24) 


From x G D N , we conclude that T <N, (x) e IK". It follows from ( |2.22[ ) that T {N) (x) G D. The inequal¬ 
ity ( |2.23| ) implies E{T iN) (x)) G J since (p N : J —> J and E(x) G J. Thus we have both T (N \x ) G D 
and E(T <n> (x)) g J. Applying Lemma 2.9 to the vector T iN) (x), we get T {N '(x) G D v which proves 
( |2.24[ ). Hence, x (0) is an initial point of T (N) . □ 

Now, we are ready to state and prove the main result in this section. In the case N = 1 and 
p - oo this result reduces to Theorem 7.3 of IflOl . 
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Theorem 2.11. Let f e K[z] be a polynomial of degree n > 2, ^ e IK" be a root-vector of f, N > 1 
and 1 < p < oo. Suppose x (0> e IK" is a vector with distinct components such that 

v F(£(jc (0) )) < 2, (2.25) 


where the function E is defined by and 'P is defined by < 12.2| ). Then f has only simple zeros 
and the Weierstrass-type iteration ( | 1 ,8| ) is well-defined and converges to f with error estimates 


\\x (k+l) - £|| < 6A (N+X)k - £|| and ||x® - f\\ < - £||, 


(2.26) 


for all k > 0, where A = c/) N (E(x {0) )), 6 = fi N (E(x (0) )) and the reed functions 0,y and (Av are defined 
by Definition 2J and Definition \2.4\ respectively. Moreover, if the inequality in ( |2.25[ ) is strict, 
then the Weierstrass-type iteration converges to ij with order of convergence N + 1. 

Proof. Since the function T is increasing on M + and R is the unique positive solution of the equa¬ 
tion 'P(t) = 2, then the initial condition ( |2.25[ ) is equivalent to E(x i{)] ) 6 J, where J = [0,7?]. Now 
the statement of Theorem |2. 11 1 follows from Theorem 3.3 of lUOl . Lemma [2. 10| and Lemma [23] 

□ 


The following result is a simplified version of Theorem 2.11 It involves only the function f 
defined by ( |2.4[ ). 

Corollary 2.12. Let f 6 K[z] be a polynomial of degree n > 2, £ 6 IK" be a root-vector off, N > 1 
and 1 < p < oo. Suppose x (0> e K" is a vector with distinct components satisfying ( |2.25[ ). Then f 
has only simple zeros and the Weierstrass-type iteration ( |1.8[ ) is well-defined and converges to L 
with error estimates 

||x ( * +1) - ^|| < A N(N+1)k ||x w - f|| and ||x® - - $\ (2.27) 

for all k > 0, where A = f{E{x^)) and the real function f is defined by ( 2.4[ ). Moreover, if the 
inequality in ( |2.25[ ) is strict, then the Weierstrass-type iteration converges to T with order of con¬ 
vergence N + 1. 

Proof. It follows from Theorem |2.11[ Lcmma [23| ii) and the inequality 0 < t/r w (0 ^ 1 which holds 
for every t e [0,7?]. □ 

We end this section with a convergence theorem under initial condition in explicit form. 

Corollary 2.13. Let f 6 IK[z] be a polynomial of degree n > 2, £ e IK" be a root-vector of f, N > 1 
and 1 < p < oo. Suppose x m> e IK" is a vector with distinct components such that 


E(x (0) ) = 


,( 0 ) 




d(x ((y >) 


< 


i(2 1/n - 1 ) 

(n - Vfi/Q + 2 


(2.28) 


where the function E is defined by ( |2.1| ) and q is the conjugate exponent of p. Then the Weierstrass- 
type iteration ( |1.8[ ) is well-defined and converges to T with order of convergence N + 1 and with 
error estimates (|2.26|) and (12.271). 


Proof. It follows from Theorem 2.1 1[ Corollary |2.12|and the lower estimate in ( |2.3| ). 


□ 
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3. Semilocal convergence theorems for the Weierstrass-type methods 


In this section, we present two semilocal convergence theorem for Weierstrass-type methods 


( |1.8[ ) with computationally verifiable initial conditions and with computationally verifiable a poste¬ 
riori error estimates. We study the convergence of the Weierstrass-type methods ( |1.8| ) with respect 
to the function of initial conditions E: T> —> M + defined by 


E f (x) = 


W f (x ) 


d{x ) 


(1 < p < oo). 


(3.1) 


3.1. First semilocal convergence theorem 

Recently Proinov IfTTI have proposed a new approach for obtaining semilocal convergence 


results for simultaneous methods via local convergence results. In particular, from Theorem 2.11 


we can obtain a convergence theorem for Weierstrass-type methods ( | 1 ,8| > under computationally 
verifiable initial conditions. In what follows q is the conjugate exponent of p. 

Theorem 3.1 (Proinov IfTTlO . Let f 6 K[z] he a polynomial of degree n > 2. Suppose there exists 
a vector x € K” with distinct components such that 


E f (x ) 


Wf(x) 


d(x) 


< 


1 


(1 + yja) 2 


(3.2) 


for some 1 < p < oo, where a = (n - \) llq . In the case n = 2 and p = oo we assume that inequality 
in ( |3.2[ ) is strict. Then f has only simple zeros and there exists a root-vector £ G K' ! off such that 


|v - £|| < a(Ef(x )) ||W/(x)|| and 


£ 


d{x) 


< 


h(E f (x)), 


where the real functions a and h are defined by 

a(t ) = 2/(1 - (a - 1)1 + V(1 - ( a ~ 1)0 2 -4 1) and hit) - ta(t). 
Moreover, if the inequality ( |3.2[ ) is strict, then the second inequality in ( |3.3[ ) is strict too. 
In the sequel, we define the real function 


(3.3) 


(3.4) 


Q(l) = W(h(t)) 


(3.5) 


where T and h are defined by ( |2.2[ ) and ( |3.4[ ), respectively. 

Now we are in a position to state the main result of this paper. 


Theorem 3.2. Let f 6 lK[z] be a polynomial of degree n > 2, N > 1 and 1 < p < oo. Suppose 
x a>> G K' ! is an initial guess with distinct components satisfying 

E f (x (0) ) < 1/(1 + xfd) 2 and Q(E f (x (0) )) < 2, (3.6) 
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where the function E r is defined by ( ]3 .1 [ ), the function Q is defined by p.5| ) and a = (n - \f !q . 
In the case n = 2 and p = oo we assume that the first inequality in ( 3.6[ ) is strict. Then f has 
only simple zeros in K and the Weierstrass-type iteration ( | 1.8 [ ) is well-defined and converges to a 
root-vector £ of f with order of convergence N + 1 and with error estimate 

llx® - £11 < a(E f (x< k >)) ||W / (x®)||, (3.7) 

for all k > 0 such that E f(x <k> ) < 1/(1 + \[af and iXE Tx®)) < 2, where the function a is defined 
by ( |X4l ). 

Proof From the first inequality in p.6| ) and Theorem 3.1, we conclude that / has only simple 
zeros and there exists a root-vector £ 6 IK" of / such that 


JO) 




d(x^) 


< 


h(E f (x m )). 


From this and the second inequality in p.6| ), taking into account that V F is increasing on [0, +oo), 
we obtain 




Vj-, 

1 

o 

7K 

^— 

\ 

d(x ((y >) 

pi 


< h(E f (x (0} ))) = Q(E f (x W) )) < 2. 


.( 0 )\ 


It follows from Theorem 2.11 that the Weierstrass-type iteration (jTT8J) is well-defined and con¬ 
verges to £ with order of convergence TV + 1. It remains to prove the estimate ( |3.7[ ). Suppose that 
for some k > 0, 

E f (x (ky ) < 1/(1 + s[ci) 2 and Q(E f (x <k) )) < 2. (3.8) 

Then it follows from the first inequality in p.8[ ) and Theorem |3.1| that there exists a root-vector 
q 6 IK" of / such that 


|x® - q\\ < a(E f (x (k) )) ||W / (x®)|| and 


Jk) 


d(x (k f 


< 


h(E f (x m )). 


(3.9) 


From the second inequality in (3.9) and the second inequality in (3.8), we get 


¥ 


x® - q 

\ 

d(x (k) ) 

pi 


< 'V{h{E f {x (ky ))) = T>(E f (fi K> )) < 2. 




By Theorem 2.11 we conclude that the iteration d 1 .8| ) converges to q. By the uniqueness of the 
limit, q = f. Hence, the error estimate ( |3.7[ ) follows from the first inequality in p.9[ ). □ 

Corollary 3.3. Let f 6 K[zJ be a polynomial of degree n >2 and 1 < p < oo. Suppose x (0) 6 IK" 
is an initial guess with distinct components such that 


Ef(x m) ) 


W f (x i0) ) 


d(x^ 0) ) 


< 


5a + 6 


(3.10) 


where a = (n - \) xlq . Then f has only simple zeros in IK and the Weierstrass-type iteration ( | 1 .8[ ) is 
well-defined and converges to a root-vector of f with order of convergence N + 1 and with error 
estimate (|3.7|). 
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Proof. The proof follows from Theorem |3.2| because the initial condition ( |3.10| ) implies ( |3.6[ ). 
For simplicity, we set R = 2/(5 a + 6) and p = 1/(1 + yfa) 2 . To prove that x (0) satisfies p.6[ ) it 
is sufficient to show that R < p and Q.(R) < 2, since Q is increasing on [0, +oo]. We prove only 
Q(R) < 2. We define the function G: M+ —> M+ (the graph of G is given in Fig. [I]) by: 

G(t) = (1 + 2g(t))e tR{, \ where g(t) = 4/(3f + 8 + V9t 2 + 8t+ 16). 



Fig. 1. Graph of the function G on M+. 


By using standard arguments of calculus, it is easy to prove that G(t ) < 2 for all t e M + . It follows 
from the definition of Q and the well-known inequality (1 + t/n) n < e r (t e M + ) that 


W) = (1 + 2 g(a)) 

which completes the proof of the corollary. 


, n— 1 


1 + < G(a) < 2 


n— 1 / 


(3.11) 

□ 


3.2. Second semilocal convergence theorem 

In this section we need another theorem of Proinov mn that enables us to transform local 
convergence theorems into semilocal ones. 

Theorem 3.4 (Proinov [[TTll ). Let f 6 K[z] he a monic polynomial of degree n > 2. Suppose the 
initial guess x 6 K" with distinct components such that 


W f (x) 

d{x) 


R(l-R) 

1 +(a- 1 )R 


(3.12) 
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for some 1 < p < oo and 0</?<l/(l + s/a), where a = (n — l) 1/<? . In the case n = 2 and p = oo 
we assume that inequality in ( |3. 1 2[ ) is strict. Then f has only simple zeros and there exists a 
root-vector £ e K" of f such that 


\x - £|| < a(Ef(x )) ||VT(.r)|| and 


d(x) 


<R, 


(3.13) 


where the real function a is defined as in ( 3.4[ ). If the inequality ( 3.12 ) is strict, then the second 
inequality in ( ]3. 1 3[ ) is strict too. 

Theorem 3.5. Let f 6 K[zJ he a polynomial of degree n > 2, N > 1 and 1 < p < oo. Suppose 
;c ,(,) 6 IK" is an initial guess with distinct components such that 


E f (x (0) ) = 


W f (x {0) ) 


d(x^) 


(V2- \){a + 2 - nif/l - 1)) 


^ n 

^ 03 _ 

p (a + 2 )(a + 2 + n(a - l)(S/2 - 1)) 


(3.14) 


where a = (n - 1) 1 lq . Then f has only simple zeros in IK and the Weierstrass-type iteration ( | 1 ,8[ ) is 
well-defined and converges to a root-vector off with order N + 1 and with error estimate 

IU W - £11 < a(E f (x <k fi) ||W / (x®)||, (3.15) 

for all k > 0 such that E f(x lk> ) < % where the function a is defined by ( |3.4[ ). 

Proof. Let R = {nififl - 1 ))/(a + 2). From the well-known inequality nifij 2 - 1) < 1, we obtain 
R < 1/(1 + s/d). On the other hand 

R(l-R) 


__ nft2- l)(a + 2 - njf/l - 1)) 

1 + (a - 1)R ~ (a + 2)(a + 2 + n(a - 1)(V2 - 1)) 


Therefore, ( |3. 14[ ) can be written in the form 

lT(a: (0) ) 


d(x w >) 


< 


R(l-R) 

1 +( fl - 1 )R 


Then it follows from Theorem |3.4| that / has only simple zeros in IK and there exists a root-vector 
£ 6 IK" of / such that 


,( 0 ) _ 


£ 


d(x^) 


<R. 


Now Corollary 2.13 implies that the Weierstrass-type iteration ( |1.8[ ) converges to £ with order of 
convergence N + 1. It remains to prove the error estimate ( |3.15| ). Suppose that for some k> 0, 

W f (x (k) ) 


Ef(x (k) ) = 


d(x (k) ) 


< 3 ? = 


R(l-R) 

1 + (a - 1 )R 


(3.16) 


Then it follows from Theorem |3.4| that there exists a root-vector q 6 IK" of / such that 


U k) -q\\<a{E f {x (k) ))\\W f {x^) 

12 


and 


y(k) 


d(x (k >) 


<R. 


(3.17) 

























































From the second inequality in ( |3. 17| ) and Corollary |2.13[ we conclude that the Weierstrass-type 
iteration ( |1.8| ) converges to r\. By the uniqueness of the limit, we get tj = f Consequently, the 
error estimate ( |3.15| ) follows from the first inequality in ( |3. 1 7| ). This ends the proof. □ 

Corollary 3.6. Let f e K[z] be a polynomial of degree n > 2, N > 1 and 1 < p < oo. Suppose 
x <(>) 6 K" is an initial guess with distinct components such that 


E f (x {0) ) 


W f (x (0) ) 


d(x (Q >) 


< 


niff '.2 - 1 ){a + 1) 
(a + 2)(2 a + 1) 


(3.18) 


where a = (/? - 1) 1/V/ . Then f has only simple zeros in IK and the Weierstrass-type iteration ( | 1 ,8[ ) is 
well-defined and converges with order N + 1 to a root-vector of f with error estimate (]3. 15[). 


Proof It follows from Theorem 3.5 and the inequality ntfjl - 1) < 1. 


□ 


4. Numerical examples 

In this section, we provide three numerical examples to show the applicability of Theorem |3.2[ 
We consider only the case p = oo since the results in other cases are similar. Let / 6 C[z] be a 
polynomial of degree n > 2 and let x <{)) e C” be an initial guess. We consider the function of initial 
conditions Ej- : T> —> M + defined by 


Ej(x) 


Wf(x) 


d(x) 


\n— 1 


Furthermore, we define the real function Q as follows 

n(t) = (1 + 2ta(t))(\ + ta{t)) n 
where the function a is defined by 

a(t) = 2/(1 -(n- 2 )t + yj{\ - (n - 2)t) 2 - At). 

It follows from Theorem |3.2| that if there exists an integer m > 0 such that 

E f (x {my ) <p = 1/(1 + V/i - 1 ) 2 and Ll{E f {x {m) )) < 2, 


(4.1) 


(4.2) 


(4.3) 


(4.4) 


then / has only simple zeros and the Weierstrass-type iteration ( |1.8[ ) starting from x (i)> is well- 
defined and converges to a root-vector ^ of / with order of convergence N + 1. Besides, the 
following a posteriori error estimate holds: 


\\x {k) - flloo < e k , where s k = a(E f (x w )) \\W f {x { 


ik)\ 




for all k > m such that 


E f (x {ky ) < p and Q(E f (x w )) < 2. 


ikh 


(4.5) 

(4.6) 


13 




































In the examples below, we apply the Weierstrass-type methods (jT78j) for some N > 1 using the 
stopping criterion 

s k < 10~ 15 (k > m) (4.7) 


together with (J4T6J). For given W > 1 we calculate the smallest m > 0 which satisfies the conver¬ 
gence condition ( |4.4| ), the smallest k > in for which the stopping criterion ( |4.7[ ) is satisfied, as well 
as the value of s k for the last k. From these data it follows that: 1) / has only simple zeros; 2) 
the /Vth Weierstrass-type iteration (| 1 ,8[) starting from x m> is well-defined and converges with order 


V + 1 to a root-vector of /; 3) at kth iteration the zeros of / are calculated with an accuracy at 
least e k . 

In Table[2]the values of iterations are given to 15 decimal places. The values of other quantities 
(jj. Ef(x (m) ), etc.) are given to 6 decimal places. 


Example 4.1. Let us consider the polynomial f(z) = z 3 - z and the initial guess 


JO) _ 


= (1.74,1.75,-3.49) 


which are taken from Hopkins et al. D21 and Niell [[63l. We have p = 0.171572. For a given N, it 
can be seen from Table |T| the value of m which guarantees that the Weierstrass-type iteration ( |1.8[ ) 
starting from x ((y> is well-defined and converges to a root-vector of / with order of convergence 
N + 1, the value of k for which the stopping criterion ( |4.7[ ) is satisfied and the value of the error 
estimate s k which is guaranteed from Theorem 3.2 For example, for N = 100 at the second 
iteration (m = 2) we have proved that the method is convergent with order of convergence 101 
and at the third iteration (k = 3) we have calculated the zeros of / with accuracy less than 10 -523 . 
Moreover, at the fourth iteration we have obtained the zeros of / with accuracy less than 10 -5290() . 
The numerical results for N = 100 are shown in Table [2] 


Table 1 


Values of m, k and e k for Example |4.l[ 


N 

m 

E f (x {m) ) 

Q(E f (x (my )) 



k 

s k 


£ k +i 

1 

12 

0.029714 

1.131702 

3.311488 x 

io- 2 

16 

5.496409 x 

10- 2 6 

3.000715 x 10" 51 

2 

6 

0.007688 

1.031545 

7.903736 x 

10~ 3 

8 

2.463566 x 

10~ 21 

7.688556 x 10~ 63 

3 

6 

0.000216 

1.000867 

2.169611 x 

io- 4 

8 

1.692612 x 

IO- 5 9 

8.138142 x IO" 236 

4 

4 

0.007479 

1.030664 

7.656408 x 

10~ 3 

6 

2.712088 x 

10-66 

1.252586 x IO" 330 

5 

6 

0.000000 

1.000000 

3.741978 x 

10~ 8 

7 

1.837441 x 

10- 45 

3.058350 x IO' 269 

6 

4 

0.000361 

1.001445 

3.613767 x 

io- 4 

5 

7.021265 x 

10- 29 

1.900890 x IO" 199 

7 

3 

0.016712 

1.070710 

1.766014 x 

io- 2 

4 

5.881957 x 

10- 17 

1.306375 x 10~ 131 

8 

4 

0.000000 

1.000000 

6.811047 x 

10~ n 

5 

1.439954 x 

10- 95 

1.144468 x IO" 857 

9 

3 

0.013852 

1.058033 

1.387643 x 

io- 2 

4 

2.122314 x 

10- 19 

1.503595 x 10" 187 

10 

4 

0.002015 

1.008114 

2.019382 x 

10~ 3 

5 

1.020330 x 

10~ 36 

2.321516 x 10‘ 402 

100 

2 

0.000006 

1.000026 

6.628377 x 

10~ 6 

3 

2.609028 x 

10- 524 

3.867338 x IO" 52901 
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Table 2 


Numerical results for Example |4.l| in the case N = 100. 


iteration 

4 ° 


x (k) 

k = 

0 

1.74 

1.75 

-3.49 

k = 

1 

1.149415748340902 

1.975676419092484 

-2.359878141616537 

k = 

2 

0.999998661360835 

-0.000006628312624 

-1.000004865683659 

k = 

3 

1.000000000000000 

0.000000000000000 

- 1.000000000000000 


Example 4.2. Consider the polynomial f(z) = z 7 - z 5 - 10z 4 - z 3 - z + 10 and the initial guess 
x (0) = (2.3 + 0.1/, 1.2 + 0.2/, -0.8 - 0.2/, 0.1 + 1.3/, -0.2 - 0.8/, -1.2 + 2.2/, -1.2 - 1.8/) 


which are taken from Sakurai and Petkovic [|T6ll . In this case jj = 0.084040. The results for this 
example are presented in Table [3] For example, we can see that for N = 4 at the first iteration 
(in = 1) we prove the convergence of the method, and at the third iteration (k = 3) we obtain the 
zeros of / with accuracy less than 10 -58 . 

Table 3 


Values of m, k and for Example |4.2[ 


N 

m 

E f { ;c (m) ) 



k 

£k 

£/c+i 

1 

2 

0.007526 

1.064790 

1.116392 x 10- 2 

5 

1.796060 x 10- 17 

2.792108 x lO" 34 

2 

1 

0.035532 

1.427605 

6.352229 x 10" 2 

4 

1.209144 x 10- 39 

7.010810 x 10- 118 

3 

1 

0.013767 

1.126494 

2.129981 x 10- 2 

3 

2.368469 x 10~ 31 

4.291912 x lO' 123 

4 

1 

0.004823 

1.040419 

6.681020 x 10- 3 

3 

1.000227 x 10- 59 

8.418384 x 10“ 297 

5 

1 

0.001903 

1.015502 

2.840694 x lO" 3 

2 

2.619223 x 10~ 17 

7.631970 x KT 101 

6 

1 

0.000695 

1.005604 

9.366066 x 10~ 4 

2 

1.157166 x 10- 22 

9.947018 x 10~ 156 

7 

1 

0.000253 

1.002031 

3.750097 x 10~ 4 

2 

1.968419 x 10" 29 

1.470015 x 10" 231 

8 

1 

0.000107 

1.000862 

1.444295 x 10“ 4 

2 

3.245945 x 10~ 36 

1.808200 x lO' 322 

9 

1 

0.000038 

1.000306 

5.655465 x 10~ 5 

2 

9.224622 x 10~ 45 

3.800354 x 10“ 443 

10 

1 

0.000015 

1.000124 

2.091765 x 10‘ 5 

2 

1.833150 x 10~ 53 

1.621635 x KT 583 

100 

1 

0.000000 

1.000000 

1.325425 x 10“ 40 

1 

1.325425 x 10“ 40 

1.089487 x lO" 4036 


Example 4.3. In 1997, Wang and Zhao [|T7l studied the convergence behavior of the Weierstrass 
method d 1,3[ ) for the polynomials 


m 


,20 


- 1 and /(z) 


,30 


- 1 . 


(4.8) 


In their paper they say: “... when we intended to solve the equation z 20 - 1 = 0, we could not find 


any suitable initial values to use algorithm ( ]1,3| )”. Actually, this is not true. There are many initial 
approximations for which the Weierstrass method ( |1.3[ ) is convergent for the polynomials ( |4.8[ ). 
We have made two types of experiments. 
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First, we compute the zeros of the polynomials (|4.8[) by using the Weierstrass algorithm (| 1 .3[) 


with 1000 initial guesses x' 


.( 0 ) 


such that |U (0, ||oo < 2 given randomly. By Theorem 3.2 (N = 1), 


we obtain that the Weierstrass method starting from each of these random initial approximations 
is well-defined and convergent to a root-vector of /. 

Second, we compute the zeros of the polynomials ( |4.8[ ) by using the Weierstrass-type iterations 
3]) for some N > 1 (including the case A = 1) with Aberth’s initial approximation x i0) e C" 


40 ); 


( 1 . 


given by (see flU) 


x[° } = r 0 exp (id v ), 6> v =^(2v-^j, v=l,...,n, (4.9) 

where n is the degree of the corresponding polynomial. We have made these experiments for 
r 0 = 1,1.1,1.2,..., 1.9,2. Again, in all cases, we prove the convergence of the methods. In this 
example, we present the results for the second type of experiments for r 0 = 2. For the polynomial 
f(z) = z 20 - 1 we have /j = 0.034821 and the obtained results can be seen in Table |4j For example, 
we can see that for N = 61 at the seven iteration we have calculated the zeros of / with accuracy 
less than 10 -14153 . 

Table 4 


Values of m, k and s k for Example [O] for f(z) = z 20 - 1. 


N 

m 

E f (x im) ) 

Q(E / (x (m) )) 


k 

Sk 

Sk+ 1 

1 

16 

0.005454 

1.135937 

1.906753 x 10- 3 

19 

5.251672 x 10“ 16 

2.620105 x 10- 30 

2 

10 

0.008641 

1.241514 

3.249990 x 10~ 3 

12 

6.054274 x 10~ 16 

2.002780 x 10~ 44 

3 

8 

0.006432 

1.165842 

2.298445 x 10~ 3 

10 

3.924632 x 10" 29 

2.034074 x 10- 111 

4 

7 

0.003429 

1.079931 

1.147442 x 10- 3 

9 

1.568679 x 10" 51 

7.736874 x 10" 251 

5 

7 

0.000000 

1.000000 

1.310563 x 10~ 8 

8 

3.920705 x 10- 43 

2.810626 x 10" 250 

6 

6 

0.000465 

1.009907 

1.469386 x 10~ 4 

7 

1.026738 x 10" 21 

8.842207 x 10" 142 

7 

6 

0.000000 

1.000006 

9.113539 x 10- 8 

7 

3.323098 x 10~ 50 

1.038511 x 10- 389 

8 

5 

0.014073 

1.494951 

6.079699 x 10~ 3 

7 

2.518063 x 10- 112 

2.700157 x 10~ 997 

9 

5 

0.001649 

1.036367 

5.324415 x 10” 4 

6 

8.150179 x 10- 25 

8.150497 x 10“ 233 

10 

5 

0.000075 

1.001583 

2.357206 x lO" 5 

6 

7.347516 x 10- 42 

2.017354 x 10“ 443 

61 

5 

0.000069 

1.001472 

2.192754 x 10- 5 

6 

5.604020 x 10- 230 

1.117175 x 10- 14154 

100 

3 

0.000000 

1.000000 

4.366726 x 10“ 17 

3 

4.366726 x 10“ 17 

2.679890 x 10- 1555 

101 

3 

0.000000 

1.000000 

1.612383 x 10‘ 17 

3 

1.612383 x 10- 17 

8.163089 x 10" 1615 


In the Figure [2| we present the trajectories of approximations generated by the methods ( |1.8[ ) 
for N = 1 after 19 iterations and N = 61 after 6 iterations. The trajectories of approximations 
generated by other Weierstrass-type methods Q1.8|) are similar either to N = 1 or to N = 61. 
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N = 1 (it = 19) 



N = 61 (jfc = 6) 



Fig. 2. Trajectories of approximations for the polynomial f(z) = z 20 - 1. 


The situation is similar for the polynomial /(z) = z 30 - 1. In this case /i = 0.024527 and the 
obtained numerical results are presented in Table [5} For example, we can see that for N = 101 at 
the sixth iteration we have obtained the zeros of / with accuracy less than 10 -95810 . 


Table 5 

Values of in, k and s k for Example |4.3| for f(z) = z 30 - 1. 


N 

m 

E f {x {m) ) 

Q(E f {x (my )) 


k 

£k 

Sk+ 1 

1 

23 

0.004903 

1.193434 

1.196341 x 10~ 3 

26 

1.664050 x 10~ 16 

4.015143 x 10“ 31 

2 

15 

0.000303 

1.009546 

6.408814 x 10- 5 

17 

3.307885 x 10" 29 

7.610048 x 10- 84 

3 

12 

0.000132 

1.004131 

2.780420 x 10~ 5 

14 

3.153464 x 10" 56 

3.014782 x 10- 219 

4 

10 

0.003933 

1.147300 

9.286229 x 10‘ 4 

12 

5.264378 x 10" 50 

1.787341 x 10- 242 

5 

9 

0.003966 

1.148751 

9.371733 x 10~ 4 

11 

4.726532 x 10~ 71 

7.146541 x 10- 417 

6 

9 

0.000000 

1.000000 

3.581766 x 10- 9 

10 

7.028904 x 10- 53 

7.878100 x 10“ 359 

7 

8 

0.000337 

1.010607 

7.117318 x 10- 5 

9 

8.244856 x 10~ 26 

2.877668 x 10~ 193 

8 

8 

0.000000 

1.000000 

4.511178 x 10- 9 

9 

1.511999 x 10- 66 

8.070636 x 10- 584 

9 

7 

0.006793 

1.299041 

1.773403 x 10~ 3 

8 

1.009302 x 10“ 18 

3.108310 x 10- 170 

10 

7 

0.000195 

1.006113 

4.110752 x 10 -5 

8 

2.188233 x 10- 37 

2.263137 x 10“ 392 

101 

4 

0.000000 

1.000000 

4.263459 x 10-“ 

5 

3.419093 x IQ" 941 

5.715983 x IQ" 95811 


Remark 4.4. Let us note that the numerical results of Examples |4.1| |4.2| and 4.3 can also be 
obtained by each of Theorem |3.5[ Corollary |3 .3 1 and Corollary |3.6| 
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